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ABSTRACT 

Context. The prospects for using asteroseismology of rapidly oscillating Ap (roAp) stars are hampered by the large uncertainty in 
fundamental stellar parameters. Results in the literature for the effective temperature (T eff ) often span a range of 1000 K. 
Aims. Our goal is to reduce systematic errors and improve the T e g calibration of Ap stars based on new interferometric measurements. 
Methods. We obtained long-baseline interferometric observations of j} CrB using the CHARA/FLUOR instrument. In order to disen- 
tangle the flux contributions of the two components of this binary star, we additionally obtained VLT/NACO adaptive optics images. 
Results. We determined limb darkened angular diameters of 0.699 ± 0.017 mas for ji CrB A (from interferometry) and 0.415 ± 
0.017 mas for ji CrB B (from surface brightness-color relations), corresponding to radii of 2.63 ± 0.09R o (3.4% uncertainty) and 
1.56 ± 0.07 R G (4.5 %). The combined bolometric flux of the A+B components was determined from satellite UV data, spectropho- 
tometry in the visible and broadband data in the infrared. The flux from the B component constitutes 16 ± 4 % of the total flux and 
was determined by fitting an ATLAS9 model atmosphere to the broad-band NACO J and K magnitudes. Combining the flux of the A 
component with its measured angular diameter, we determine the effective temperature (A) = 7980 ± 180 K (2.3 %). 
Conclusions. Our new interferometric and imaging data enable a nearly model-independent determination of the effective temperature 
of fi CrB A. Including our recent study of a Cir, we now have direct T cff measurements of two of the brightest roAp stars, providing a 
strong benchmark for an improved calibration of the r cff scale for Ap stars. This will support the use of potentially strong constraints 
imposed by asteroseismic studies of roAp stars. 

Key words, techniques: interferometric, stars: chemically peculiar, stars: fundamental parameters, stars: individual: /? CrB, a Cir, 
y Equ, 10 Aql 



^ ■ 1. Introduction 1.1. A brief history of f3 CrB 
5-j | 

Photometric and spectroscopic determinations of the effective p CrB (HD 137909, HIP 75695) is one of the brightest, coolest, 
temperatures of Ap stars are affected by systematic errors. This and best-studied magnetic Ap stars. The literature on the star 
has been corroborated by asteroseismic data of rapidly oscillat- is extensive and we only mention a few of the most impor- 
ing Ap (roAp) stars in general and, more recently, by the first in- tant res ults here. P CrB has been c lassified as type A9 Sr Eu Cr 
terferomet ric determination o f the angular diameter of the roAp star by |R e nson & Manfroldl (120091). Its binary nature was first 
star a Cir dBruntt et alj|2008|). Unfortunately, the intriguing as- suggestedbv TCampbell & Moord (119071) . and recent determina- 
teroseismic potential offered by roAp stars is strongly compro- tions of its orbital elem ents were obtained bv lTokovininl (fl98l) 
mised by the presence of these systematic errors. We therefore and lNorth et all (fT^l . From sp eckle interferometri c measure- 
seek to make direct measurements of the radii and effective tern- ments using narrow-band filters, iHorch et alJ d2004h measured 
peratures of a number of Ap stars using interferometric and spec- the magn itude difference to be 2.37 mag at 55 1 nm and 1 .99 mag 
trophotometric measurements. We will first give a brief sum- at 503 nm. In the analysis p resented in Sect. [3] we retain the 
mary of the properties of /3 CrB before describing our obser- orbital elements obtained by iTokovininl (1 19841) . as they are in 
vations, data reduction (Sect. H and analysis (Sect. [3]). significantly bette r agreement with our NACO astrometry than 

Nort hetaiT (fT998). 



Send offprint requests to: H. Bruntt iNeubaueri dl944 suggested that a third body could be 

Based on observations made with ESO telescopes at the La Silla present in the system, causing r adial velocity variation s 



Paranal Observatory, under ESO DDT program 281.D-5020(A). with a period of » 32 1 days, but [Oetken & Orwertl dl984l) . 

Correspondence to: bruntt@phys.au.dk iKamper et alJ d 1990b and Soderhielm ( 1999j) excluded this pos- 
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sibility. Recently. iMuterspaug h et al. (2006) established an up- 
per limit of * 10 to 100 My (depending on the orbital period) 
for a possible substellar tertiary fr om differential interferomet- 
nc astrometry. iTrilling et all (2007) searched for 24 and 70 /im 
infrared excess around p CrB using Spitzer but did not find any. 
Interestingly, the Spitzer flux they obtained is significantly below 
the expected flux at 24 p.m, and slightly lower (although compat- 
ible) at 70 /mi. This result could be due to the chosen physical 
parameters for their stellar atmosphere model. In the following, 
we will therefore consider that /3 CrB is a binary system. 

1.2. Is p CrB a pulsating star? 

Early photometric searches for pu lsation in /J CrB (e.g. 
Weis s & Schneiderll989t iKreidll 1 99 lh gave null results and this 
contributed to the discussion of the existence of non-oscillating 
Ap stars ("noAp";|Kurtl[l989). This has changed since the ad- 
vent of large telescopes and ultra-stable spectrog raphs. Based on 
spectr oscopic time series of a single Fe line, iKochukhov et al.l 
(2002) claimed the first possible detection of a pulsation mode 
in P CrB with a period o f 11.5 min. This result was questioned 
by Hatzes & Mkrtichianj d2004l) and was also not confirmed 
by iKurtz et al.l d2007l) . However, the goo d agreement between 
the independent spe ctroscopic stu dies of [H atzes & Mkrti chianl 
(12004 . IKurtz et al.1 (120071) and IKochukhov eta l. (2008) con- 
firmed that p CrB is indeed an roAp star with a single known 
low-amplitude mod e with period 16 .2 min. The most robust re- 
sult was found bv lKurtz et alj d2007l) . who used 2 hours of high- 
cadence time-series spectra obtained with VLT/UVES. They 
detected a single oscillation frequency at 1 .03 1 mHz (P - 
16.2 min) with an amplitude of 23.5 ± 2.4kms~' in the Ha line 
and a higher amplitude in the cesium lines. Unlike most roAp 
stars, variation was observed only in singly-ionized rare-earth 
elements, b ut not doub l y ioni zed lines. The abundance analy- 
sis done bv lKurtz et alj d2007l) on their averaged spectru m con- 
firmed earlier investigations by Rvabchikov a et al] d2004l) . These 
analyses show that p CrB has an overabundance of rare-earth el- 
ements but only by about 1 dex, contrary to the 2-3 dex seen in 
most roAp stars. 

2. Observations and data reduction 

2. 1 . VLT/NACO adaptive optics imaging 

We observed p CrB on 16 Jun e 2008 using the N asmyth 
Adaptive Optics System (NAOS; iRousset et all 120031) of the 
Very Lar ge Telescope (VLT) , coupled to the CONICA infrared 
camera dLenzen et al.1 LL998), abbreviated as NACO. We se- 
lected the smallest ava ilable pixel scale of 13.26 + 0.03 mas/pix 
dMasciadri et al.ll2003l) . giving a field of view of 13.6"xl3.6". 
This small scale resulted in good sampling of the Point Spread 
Function (PSF). We employed the J and K filters of NACO, with 
respective bandpasses of 1.265+0.25 /mi and 2.18+0.35 /mi, to- 
gether with a neutral density filter (labeled "ND2_short", trans- 
mission m 1.5 %) to avoid saturation of the detector. The trans- 
mission curves of these filters are available on the NACO instru- 
ment web pag4B 

We obtained 20 images in the J band, and 40 images in the 
K band, each with an exposure time of 0.35 s. This is the min- 
imum full-frame integration time of CONICA. The / band im- 
ages were collected during « 2 minutes around UT01:53:31, and 
the K images in « 4 minutes around UT01:58:20. During these 

1 http://www.eso.org/instruments/naco/inst/filters.html 
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Fig. 1. Average NACO image of p CrB A and B in th e K band, 
together with the binary orbit from iTokovininl d 1984b . The po- 
sitions of p CrB A and B measured on our NACO image are 
marked with a "★" and "+" symbol, respectively. The scale of 
the two axes are in arc seconds, relative to p CrB A. The linear 
scale in AU is based on the Hipparcos parallax. 



observations the DIMM seeing at Paranal in the visible was good 
(0.6 - 0.7"), resulting in a high Strehl ratio (» 45 - 50%). The 
raw images were dark-subtracted, flat-fielded (using lamp flats), 
and corrected for bad pixels using IRAlQ On the NACO images, 
we measured both the differential photometry and the differen- 
tial astrometry of p CrB B relatively to p CrB A taken as the 
reference. 

To measure the relative astrometry, we treated each image 
separately using the Yoriclfl software package. We used a clas- 
sical x 1 minimization to fit an extracted sub-image of p CrB A 
(with a size of 9 x 9 pixels) at the position of the fainter com- 
ponent B. The interpolation of the shifted image of A was done 
in Fourier space. The adjusted parameters were the relative posi- 
tions c/(RA) and af(Dec), the flux ratio and the background level, 
although we used only the relative separations for our astro- 
metric analysis. In order to estimate the associa ted error bars, 
we used the bootstrapping technique described bv lKervella et al.l 
(2004a). This technique is also called "sampling with replace- 
ment" and consists of constructing a hypothetical, large pop- 
ulation derived from the original measurements and estimate 
the statistical properties from this population. The technique 
allows us to compute meaningful confidence intervals without 
any assumption on the properties of the underlying population 
(e.g. a Gaussian distribution). We validated the adopted Fourier 
interpolation method by comparing the results with a simple 
Gaussian fit of the two PSF cores. The two methods yield ex- 
actly the same relative positions (within 150 //as), although the 
Gaussian fit has a slightly larger dispersion due to the mismatch 
of the slightly seeing-distorted PSF and the Gaussian function. 

2 IRAF is distributed by the NOAO, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under co- 
operative agreement with the National Science Foundation. 

3 http://yorick.sourceforge.net/ 
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Table 1. Interferometric calibrators selected from M erand et alj 
(2005). 





Spectral 


K 


UD diameter 


HD 


type 


[mag] 


in K band [mas] 


147266 


G8III 


3.8 


0.785 ±0.011 


145457 


KOIII 


4.1 


0.693 ± 0.009 


108123 


KOIII 


3.7 


0.929 ±0.012 



We obtained the following vector separations along the RA and 
Dec directions of B relatively to A, for the epoch of the observa- 
tions (MJD 54633.08): 

</(RA) = 177.84 + 0.09 + 0.40 mas, (1) 

rf(Dec) = -208.38 + 0.06 + 0.47 mas. (2) 

The two stated error bars are the statistical and systematic un- 
certainties, respectively. The latter includes the pixel scale un- 
certainty and the detector orientation uncertainty. These values 
correspond to a separation r and position angle a (east of north) 
of: 

r = 273.95 + 0.63 mas, (3) 

a = 139.52 ± 0.13 degrees. (4) 

These measurements were done on the K images since they have 
the best Strehl ratio. In the J band, the Strehl ratio was lower and 
more unstable, resulting in a significantly variable background 
from A to B. Although its average value is not a concern, its 
slope tends to slightly shift the average apparent position of B, 
by -0.2 and -1.4 mas in RA and Dec, respectively (towards the 
lower left quadrant of Fig.Q~|i. The average NACO ima ge in the K 
band i s presented in Fig.Q] together with the orbit by Tokovinin 
(119841) . Our astrometric measurement falls on the predicted orbit 
within only 7 mas. Note that th e reference epoch of the orbital 
elements by iNorth et al] d!998l) appears to be late by approxi- 
mately 300 days. 

The photometry was obtained in two steps: (1) we obtained 
the combined photometry of the two stars; (2) we then measured 
the differential flux of B relative to A. We will discuss these steps 
in the following. 

(1) The combined ADU count was computed from the 
NACO images by using a large aperture enclosing the full PSFs 
of the two stars. It was then converted to magnitude using the 
Zero Points obtained routinely by the observatory on the same 
night: ZP(J) = 24.062 + 0.048 and ZP(K) = 22.985 + 0.078, 
and an attenuation of 4.81 ± 0.03 mag for the neutral density 
filter. These zero points are not corrected for atmospheric ab- 
sorption, but as they were obtained at low airmass (« 1.15), we 
neglect the atmospheric absorption of m 0.01 mag in J and K. 
We co rrected the atmospheri c absorption using the standard val- 
ues bv lNikolaev et alj d2000h . namely 0.092 mag/AM (relative to 
unit airmass) for J and 0.065 mag/AM for K, for our observation 
airmass of 1.71. We obtain: 

m y (A + B) = 3.28 + 0.07, (5) 
m K (A + B) = 3.22 ± 0.08. (6) 

(2) The differential photometry was obtained slightly differ- 
ently, since the diffuse background of B CrB A tends to contam- 
inate the flux of star B (but the reverse effect is negligible). We 
first computed aperture photometry of A on the average J and K 



Table 2. Journal of observations. 





MJD 


B 


PA 


V 1 % 


V-(*) % 




-54000 


(m) 


(deg) 


observed 


corrected 


A 


604.388 


324.09 


6.17 


33.2 ±1.6 


54.7± 2.7 


A 


604.408 


323.61 


1.40 


31.5 ±1.4 


52.0± 2.3 


A 


604.433 


323.88 


-4.70 


29.0 ±1.2 


47.8± 2.0 


A 


604.488 


327.21 


-17.40 


34.1 ±3.6 


56.3± 5.9 


A 


605.296 


245.70 


-39.57 


45.9 ±3.1 


75.8± 5.1 


A 


605.326 


248.92 


-45.62 


44.7 ±2.2 


73. 8± 3.7 


B 


604.398 


323.76 


3.65 


3.44±0.25 


70.0± 5.2 


B 


605.254 


238.03 


-29.58 


4.76±1.38 


97.1±28.0 



Notes: B is the projected length of the interferometric baseline, 
PA is the projection angle of the baseline, V 2 is the observed 
squared visibility and V 2 (*) is the squared visibility corrected 
for the presence of the other component (see Eq.ll2b. 



images using very small aperture radii of 3 pixels in the J band 
and 4 pixels in the K band (contamination is lower in K). We 
calculated the median background value in concentric rings cen- 
tered on A. This value was then subtracted from component B's 
flux. This allowed us to subtract the diffuse light from the PSF 
wings of A at the position of B. We checked that the residual 
background around B was negligible. We then integrated the flux 
of B on the ring-median-subtracted image using the same aper- 
ture radius as for A. We obtain the following flux ratios of each 
star relative to the total of the two, i.e. p(*) = f(*)/f(A + B): 

pj (A) = 0.7957 + 0.0084, Pj (B) = 0.2043 ± 0.0021, (7) 

p K (A) = 0.7785 + 0.0002, p K (B) = 0.2215 + 0.0002. (8) 

Note that the quoted uncertainties are statistical errors and do not 
include possible flat-fielding errors. From the combined magni- 
tudes determined above, we obtain the individual magnitudes of 
B CrB A and B: 

mj (A) = 3.54 + 0.07, mj (B) = 5.00 + 0.07, (9) 

m K (A) = 3.50 ± 0.08, m K (B) = 4.86 + 0.08. (10) 

The individual 7, K magnitudes have large uncertainties, but we 
stress that we only use the values of p for the interpretation of 
our interferometric data (Sect. 12.2) . and they are known with a 
much higher accuracy. 

2.2. CHARA/FLUOR interferometry 

Our interferometric observations of B CrB took place on 17- 
18 May 2008 in the near i nfrared K' band (1.9 < A < 
2.3 jum) at the CHARA Array dten Brummelaar et al.ll2005l) us- 
ing FLUOR (the Fiber Linked Unit for Optical Recombination; 
ICoude du Foresto et alJ 120031) . We used the FLUOR Data 
Reduction Software (DRS: ICoude Du Foresto et al.1 119971; 
iKervella et al.l l2004a: M erand et al.ll2006l) to extract the squared 
instrumental visibility of the interference fringes. We used three 
different interferometric calibrators in order to calibrate the visi- 
bilities on sky. Their properties are listed in TableQ] We note that 
the angular diameters of the calibrator stars are comparable to or 
larger than the target star. Therefore they contribute significantly 
to the uncertainty of the angular diameter measurement. The cor- 
responding systematic uncertainties were propagated into the fi- 
nal angular diameter uncertainties. 

The light from both stars of B CrB is injected simultane- 
ously in the FLUOR fibers, since the acceptance angle is 0.8" 
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on the sky. However, due to their on-sky separation of * 0.3" 
(cf. Eq. 0, two fringe packets are formed at different optical 
path differences. 

For this reason, we have to correct our measured visibility 
for this effect. The monochromatic visibility of the binary is: 

V = PK (A) V(A) + PK (B) V(B)e 2inBT/A , (11) 

where Pk(A) and Pk(B) are the relative fluxes of A and B, V(A) 
and V(B) are the individual visibilities, B the baseline vector, 
r the angular separation between A and B, and A is the wave- 
length. Because we observed over a relatively broad wavelength 
range with FLUOR, and since the binary is well-resolved by our 
baselines, the fringes appeared as two distinct fringe packets. 
Moreover, FLUOR has a limited window in terms of optical path 
difference, corresponding to a limited field of view. For (3 CrB 
A+B, the two fringe packets are not present in a single fringe 
scan. Hence, the squared visibility measured by FLUOR in the 
case of A is 

V 2 = p K (A) 2 V 2 (A). (12) 

In the case of the observations of B, this multiplicative factor 
Pk(&) 2 is small (due to the faintness of the star) and causes an 
amplification of the error bars on the true visibility V 2 (B) for 
our second observation of this star. After reducing and calibrat- 
ing the data with the DRS pipeline for each component, we have 
corrected the visibilities for this effect and then derived the an - 
gular diameters using limb-darkened models from lClarel (|2000). 
This leads to the following angular diameters: 

LD (A) = 0.699 ± 0.017 mas (2.4 %), (13) 

0ld(B) = 0.515 + 0.054 mas (10.5%). (14) 

The angular diameter of fi CrB B is significantly more uncertain 
than that of A, due to the very low apparent visibility ("V 2 ob- 
served" in Tabled of its interference fringes "on top" of the in- 
coherent flux from A. The large multiplicative factor p/r(B) ss 20 
(Eq. [T2l also applies to the associated V 2 error bar. Only two 
measurements of B could be derived from our data, and one of 
these points (MJD 54604.398) dominates the fit. In the case of A, 
our six data points are in fair agreement with the best-fit model 
with a reduced^ 2 of 2.0. 

2.3. Surface brightness-colour predictions 

We can compare the measured angular diameters of ft CrB A and 
B with the predictions from the su rface brightnes s -colour (here- 
after SBC) relations calibrated by iKervella et al.l (l2004bl) using 
their (V, V — K) relation. The K band magnitudes were obtained 
in Sect. 12.11 We derive the V band magnitudes from the total 
magn itude of the system of my(A + B) = 3.67 + 0.05 dRufenerl 
1988), and the magnitude difference Am = 1.99 + 0.10 measured 
bv lHorch et al.l (120041) by speckle interferometry at 503 nm. We 
have adopted estimated uncertainties on my and Am since they 
are not given explicitly in the references. This gives the compo- 



nent magnitudes in V: 

m v (A) = 3.83 + 0.05, m v (B) = 5.82 + 0.10, (15) 
and the predicted photospheric angular diameters: 
0ld(A) = 0.696 + 0.028 mas (4.0 %), (16) 
0ld(B) = 0.415 + 0.017 mas (4.1 %). (17) 



The predicted and measured angular diameters of j3 CrB A are 
therefore in excellent agreement, while they are in satisfactory 
agreement for B (difference of 1.7 cr). We note than increas- 
ing the uncertainty on my by a factor two does not significantly 
change the uncertainty of the angular diameters. 

2.4. Linear photospheric radii 

The original Hipparcos p arallax of ft CrB is n = 28.60+0.69mas 
dPerrvman & ESAll 19971) . consistent with the new reduction by 
Ivan Leeuwenl d2007l) of 7T = 29.17 + 0.76 mas. However, as the 
new reduction is not corrected for binarity effects, we adopt the 
original Hipparcos parallax. 

For /3 CrB A, the angular diameter measurement presented 
in Sect. 12.21 represents a significant improvement in accuracy, 
by a factor of 1 .6, over the surface brightness-colour estimate 
of Sect. 12.31 For the B component, this is not the case, as the 
visibility measurement is made particularly difficult by the pres- 
ence of the brighter A component. For the subsequent analysis 
presented in Sect. [3] we therefore choose to adopt for j3 CrB A 
our direct interferometric angular diameter measurement, while 
for B we will use the SBC estimate computed from our K band 
NACO photometry. This gives the following linear radii: 

flinterf(A) = 2.63 ± 0.09 R Q (3.4 %), (18) 
/?sbc(B) = 1.56 ± 0.07 R (4.5 %). (19) 



3. The effective temperatures and masses of fi CrB 

In the following we will determine the effective temperatures 
and luminosities of the components of /3 CrB using two methods. 
The first method (Sect. 13.11 ) relies on the bolometric correction 
(model-dependent), while the second method (Sect. 13.21 ) is only 
weakly model-dependent. We will then compare the radius and 
J e ff of the components with a grid of isochrones to determine 
their approximate age and evolutionary masses (Sect. [331 ). 

3.1. Luminosity and effective temperature from angular 
diameter + BC + parallax 

We us ed the bolometric corrections (BCy) from iBessell et al.l 
d!998l) . For the measured V - K values of 0.36 + 0.09 and 
0.99 + 0.09, for A and B we getBCy(A) = 0.04, BCy(B) = 0.01, 
BCk(A) = 0.39, and BQc(B) = 1.01 (with logg = 4). We as- 
sume the uncertainties on the BCs are 0.02 mag. We therefore 
obtain from the V band photometry: 

Wb i(A) = 3.87 + 0.05, m bol (B) = 5.83 + 0.10. (20) 

The same computation with the K band magnitudes gives iden- 
tical values within 0.01 mag. We thus obtain the following 
bolometric luminosities, assuming the Hipparcos parallax and 
M bo i = 4.75 for the Sun (recommendation by IAU 1999): 

L BC (A) = 27.6 +1.8 L , L BC (B) = 4.5 + 0.5 L . (21) 

We can now use the radii determined in Sect. 12.41 to de- 
rive the effective temperatures of the two stars through L/L Q = 
(R/Rq) 2 (Teff/Teffc) 4 , where we use the solar value T e ff ; © = 
5777K dCoxll2000l) : 

r B ff c (A) = 8 1 60 + 200 K, r e B ft c (B) = 6750 ± 230 K. (22) 
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Fig. 2. The sol id black line shows the flux d istribution for P CrB A +B, obtained by combining IUE satellite U V data, spectropho- 
tometry from lBregerl ( 1 19761) (triangles) and lAlekseeva et al] ([1996) (thick curve), and broad band fluxes from Morel & Mag nenaH 
d!978l) (box symbols). The curve peaking at 0.02 x 10 9 erg/s/cm 2 /A is the ATLAS9 model spectrum fitted to the B component, 
scaled by fitting the V and NACO J,K magnitudes (filled circles). The open circles indicate the assumed 3-sigma error bars on the 
combined flux of the components at six different wavelengths. The inset shows the details from 2000-8000 A on a linear flux scale. 



3.2. Effective temperature from angular diameter + flux 

The above method for the determination of T e g has the caveat 
that it relies on the bolometric correcti on being valid for these 
stars. The BCs from Bess ell et al. (1998) are calculated from at- 
mospheric models and do not depend on the metallicity. As a 
check, we will determine the T e ff of the A component by a di- 
rect method, meaning it will only be weakly dependent on the 
assumed model atmosphere. This is done by calculating the in- 
tegrated bolometric flux of A + B and subtracting the flux from 
the B component using an ATLAS9 model with the T e ff deter- 
mined above. 

The bolometric flux of the combined system, B CrB A+B, 
was obtained by combining data from the literature from the UV 
to the near IR as shown in Fig. [2] In the UV range we used five 
spectrograms from the Sky Survey Telescope dJamar et al.ll97 
obtained at the IUE "Newly Extracted Spectra" data archivi 
covering the wavelength interval 1150 A < A < 3350 A. We 
computed the weighted average of the spectrograms and some 
spurious data from 1150-1250 A were remo ved. In the op- 
tical range we used spectrophotometry from lAlekseeva et al.l 
(1996) and iBregerl (119761) . which cover the range 3200 A< A < 
10800 A. In addition, for the near-I R wavelengths we used th e 
broadband magnitudes VJKL from lMorel & Magnenatl (Il978l) . 
We interpolated the points between the broadband data and made 
a linear extrapolation at the end points (UV and IR ranges), al- 
though in practice the contribution is negligible. Finally, we cal- 
culated the weighted average flux, which is shown as the solid 



black line in Fig. [2] For the relative flux uncertainties we as- 
sumed 15 % in the UV, 6% in the optical, and 10% in the near 
IR. For the ranges where extrapolations were made, we doubled 
these errors. The adopted uncertainties are larger than the origi- 
nally published values. We adjusted them based on the disagree- 
ment between different sources of data in th e same wavele ngth 
ranges, i.e. the spectrop hotometric data from lBregerl (1 19761) and 
lAlekseeva etaTI d 19961) in Fig. g] 

Since B CrB is a binary system, extra care must be taken 
when computing the bolometric flux of the primary star. The 
binary has a maximum angular separation of 0.3" and all avail- 
able flux data contain the combined light of the two components. 
Since our main interest is the A component, we have to estimate 
and subtract the flux of the B component. To accomplish this, we 
fitted Kurucz models to the m j, and my magnitudes of the B 
component, taking into account the statistical errors on the mag- 
nitudes (solid circles in Fig. O. For the mj and niK magnitudes 
we adopted the measurements from NACO, while the my magni- 
tude used was that derived in Eq.[l5] The spectra for the Kurucz 
models were obtained using the IDL routine kurgetl (ATLAS9 
models) and the corresponding database of models available in 
the IUE reduction and data analysis package IUEDAC0. We 
started by converting t he my magnitude of B CrB B into flux us- 
ing the relation (6) of iRufener & Nicole j ([l988) and used this 
result to calibrate the models. We then converted the Kurucz 
fluxes at the NACO J and K central wavelengths (12,650 A and 
21,800 A) into magnitudes using m = — 2.5 xlog 10 (///o). Here m 
is the magnitude in a given filter, / is the flux at the central wave- 



4 http://sdc.laeff.inta.es/cgi-ines/IUEdbsMY/ 



5 http://archive.stsci.edu/iue/iuedac.html 
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length of that filter and /o is the standard zeroth-magnitude flux 
for the same filter. The values of /o were computed by integrating 
the flux of Vega through each of the J and K filters of the NACO 
instrument, and assuming that Vega has zero magnitude in all 
bandsS We then generated 100 values for m j and m# of/? CrB B 
by adding random fluctuations consistent with the uncertainties. 
For each set we determined the Kurucz model that best fitted 
each set of magnitudes. The average integrated flux of the 100 
fitted Kurucz models was found to be fs = (1.2 + 0.3) 10~ 7 
erg/s/cm 2 . We adopted a rather large uncertainty (25 %) since 
we only have three broadband flux measurements of the B com- 
ponent. 

The observed flux from the combined system was computed 
by integrating the black curve shown in Fig. [2] from which we 
obtained /a + fs = (7.8 + 0.4) 10~ 7 erg/s/cm 2 and the bolometric 
flux of the primary component is thus found to be /a = 6.6 + 
0.5 10~ 7 erg/s/cm 2 . We then obtained the effective temperature 
using the relation, 



Table 3. Measured quantities and derived fundamental parame- 
ters. 



—4 4 /i, i 



(23) 



where /b Q i is the bolometric flux and 6 LD is the limb-darkened 
angular diameter. Inserting the values for /a and the angular di- 
ameter from Sec. 12. 21 we obtain an effective temperature of 



(24) 



T 9 eS (A) = 7980+ 180 K. 

Combining this with the radius we get the luminosity 
L e (A) = 25.3±2.9L . 



(25) 



These values are in agreement with those in Sect. 13.11 where 
we used the (model dependent) bolometric correction. Since 
the calculation using Eq.|23]is nearly model-independent (limb- 
darkening coefficients depend on atmosphere models), we adopt 
these values as our final estimates of T e ff and the luminosity. We 
note that if we neglect the contribution to the flux from the B 
component, (A) becomes 350 K higher. 

Several determinati ons of T e g are found in the lite rature and 
we mention a few here. lKochu khov & Bagnulo (2006) used pho- 
tometric indices to determine T e g and found 7430+200 K, w hich 
is significantly lower than our value. iNetopil et alj (2008) have 
determined the T e ff of £1 CrB from three photometric systems 
(Stromgren, Geneva and Johnson) and compared these with val- 
ues in the literature. The mean value for the photometric in- 
dices is 7710 + 260 K and the mean of the literature values is 
8340 + 360 K. This is a typical example of the large scatter found 
for chemically peculiar A stars. Ho wever, since the rms scatter 
is large, the results summarized by Netop iTet al.l (120081) are in 
acceptable agreement with our new determination. It is worth 
stressing that our determination is the first that is not affected 
by photometric calibration errors or interstellar reddening, and 
is only weakly dependent on the adopted limb darkening. 

3.3. The evolutionary status and mass of fi CrB 

To investigate the evolutionary status of the two components 
of (3 CrB we have compared the ob served r e ff and radius with 
isochrones from the BASTI grid dPietrinferni et alj 120041) as 
shown in Fig. [3] To transform the mass fraction of heavy el- 
ements (Z) of the isochrones to spect roscopic IFe/H] v alues 
we used the solar value Z = 0.0156 dCaffau et alj|2009t) . i.e. 



6 We note that lBohlin & GillilanJ d2004h recently found V = 0.026 : 
0.008 for Vega. 



Parameter 


A 


B 


Ti [mas] 


ZS.OU+U.oy 


Same as A 


6 LD [mas] 


0.699±0.017 


0.415+0.017 


/ bol [10- 7 erg/s/cm 2 ] 


6.6±0.5 


1.2+0.3 


fl/Ro 


2.63±0.09 


1.56+0.07 


LIU 


25.3+2.9 


4.5+0.5 


Teff [K] 


7980+180 


6750+230 


Spectral type 


A5 


F2 


M/M Q 


2.09+0.15 


1.40+0.10 


Age [Gyr] 


0.53+0.10 


Same as A 



3.0 



2.5 




K 

a; 



2.0 



1.5 




3.94 



3.90 



3.86 



3.82 



log T e , 



Fig. 3. Radius-r e ff diagram showing the location of the two 
components of fi CrB. The thin lines are BASTI isochrones for 
[Fe/H] = +0. 10 for ages 0.5, 0.6 and 0.7 Gyr. The thick lines are 
for the sames ages but for higher metallicity [Fe/H] = +0.28. On 
some of the isochrones the box symbols mark the masses at 1.3- 
1.5 M and 2.0-2.2 M in steps of 0.1 M Q . 



[Fe/H] a; log 10 (Zbasti/Z ). We assumed an uncertainty of 
+0.002 on Z , which corresponds to +0.05 dex on [Fe/H]. 

In Fig. [3] we show two sets of isochrones with [Fe/H] = 
+0.10 and +0.28. The higher metallicity appears t o be in bet- 
ter agr eement with the location of the B component. iKurtz et al.l 
d2007b found [Fe-n/H] = +0.41 ± 0.22 from 11 lines of singly 
ionized Fe. This is the metallicity in the photosphere but we as- 
sume it represents the entire star. Since radial stratification of Fe 
is well-known to be present in roAp stars, our assumed metal- 
licity is an approximation, but it seems to be supported by the 
agreement with the location of the stars in the radius— T^g dia- 
gram in Fig. [3] With this assumption we determine the age to 
be 0.53 + 0.10 Gyr and the masses of the components to be 
Ma /M Q = 2.09 ± 0.15 and M B /M = 1.40 ± 0.10. These "evo- 
lutionary masses" ar e in good agreem ent with the dynamical 
masses determined by North^taLRl998j): M A /M = 1.87+0.13 
andM B /M = 1.356 + 0.073. 
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4. Discussion and conclusion 

We have determined the effective temperature of the primary 
component of the binary (3 CrB using a technique that is only 
weakly model-dependent. We used interferometric data to mea- 
sure the angular diameter and the fluxes were constrained using 
NACO J, K measurements of each individual component in the 
binary. 

We determined the primary component of f3 CrB to have 
T e s = 7980+1 80 K. In comparison, literature values for the com- 
bined star span 7230-8700K (considering l-cr uncertainties). 
From a similar analysis of flux data and interferom etric data the 
roAp star a Cir, we found r eff = 7420 + 170 K dBruntt et al.l 
2008). For that star the literature values also span a large range 
from 7470-8730 K. It is interesting that for /? CrB A, our re- 
sult is in the middle of the range of previous estimates while for 
a Cir the estimate is at the low end of the range. If we compare 
solely with T e g estima t es fro m the same photometric system, 
iKochukhov & Bagnulold2006l) found 7430+200 K for CrB and 
7670 + 200 K for a Cir. We must remember that the photometric 
index of ft CrB includes both components, and will always give 
a systematically low te mperature. Taking this into ac count, the 
photometric values from Kochukh ov & Bagn ulo (2006) seem to 
agree with our fundamental (i.e. model-independent) r e (j values 
within about +300 K. It will be necessary to make interferomet- 
ric measurements of several more of the brightest Ap stars to be 
able to improve the T e ff scale of these peculiar stars. 

Accurate determinations of T e g have important impact for 
the asteroseismic modelling in future work, a Cir was observed 
for 84 days with the 52 mm star trac ker on the n o w de funct 
WIRE satellite. From the light curves iBruntt et all (120091) de- 
tected five frequencies of which two had not been observed be- 
fore. These two lie symmetrically around t he well-known dom- 
inant mode at 2442 //Hz to form a triplet. iBruntt et ail {2009) 
interpreted the equidistant separation as half the large separa- 
tion. Combining this with the new T e g, the properties of the 
star could be constrained based on preliminary theoretical mod- 
elling of the observed pulsation modes. To obtain similar re- 
sults for ft CrB would be worthwhile now that all ingredients 
for the modelling of the star are available, making it the sec- 
ond roAp star with well-established fundamental atmospheric 
param eters. This would r equire an ambitious asteroseismic cam- 
paign (iKurtz et al.l l2007) using a network of telescopes with sta- 
ble spectr ographs like the Ste llar Observations Network Group 
(SONG: lGrundahl et alJl200l) . 

Our understanding of roAp stars would benefit from obtain- 
ing interferometric angular diameters of more targets. However, 
even with the most sensitive beam combiners currently available, 
only a handful are bright enough to yield a radius measurement 
to better than 2 %. Having now measured a Cir and ft CrB, we 
next propose to observe y Equ and 10 Aql. 
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